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We recall the generalization of the Feynman-Metropolis- Teller approximation for a 
compressed atom using a relativistic Fermi-Thomas model. These results within a 
Wigner-Seitz approximation lead to a new equation of state for white dwarfs and to 
a new value of their critical mass, smaller than the one obtained by Chandrasekhar. 
The possible observations of these effects in binary neutron stars are outlined. 

1 Introduction 

The problem of critical mass against gravitational collapse is one of the most 
important in relativistic astrophysics and has been the subject of a very large 
number of theoretical investigations with many observational consequences, 
the most important one being the paradigm for the identification of black holes 
in binary X-ray sources. The concept of critical mass appears in a large vari- 
ety of astrophysical settings: white dwarfs, neutron stars, galactic nuclei. The 
occurrence of gravitational collapse in galactic nuclei is still matter of great 
debate both from the observational and theoretical point of view: the physi- 
cal conditions leading to the onset of gravitational collapse, even the matter 
constituents of such a process, are still far from settled although the exis- 
tence of extremely massive relativistic systems in the range of 10 6 — 10 8 M© in 
galactic nuclei is certain. The occurrence of gravitational collapse in neutron 
stars has obtained grandiose observational verification in binary X-ray sources 
although, from a theoretical point, it still presents problems of unsurpassed 
difficulty due to theoretical and observational unknowns in the behavior of 
matter compressed at supra- nuclear densities. In this case only a firm up- 
per limit on the numerical value of the critical mass has been obtained from 
first principles: the validity of general relativity, the causality non-violation 
principle and the behavior of matter at a fiducial nuclear densityEj 

In this respect the case of white dwarfs appears most interesting: it is 
amenable to a complete theoretical treatment which is also testable by refined 
observational measurements thanks to the superb observations of binary pul- 
sars. The critical mass for white dwarfs occurs at densities relatively low in 
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comparison with nuclear densities and many of the intricacies of the behavior 
of matter in bulk at supra-nuclear densities can be completely neglected. On 
the other hand the theoretical analysis of the white dwarf critical mass can 
become the arena for a detailed treatment of gravitational and electromag- 
netic effects at relativistic regimes in the fully tested region of validity of these 
field theories: general relativity, quantum theory and electromagnetism. I will 
devote this contribution to this topic which has reached complete maturity 
and in principle can give important observational clues for the understanding 
of astrophysical high-precision measurements. 

2 Historical Background 

That white dwarfs cannot have an arbitrary large mass and that a critical 
mass against gravitational collapse must exist in their configuration of equi- 
librium reachable by increasing their central densities was clearly understood 
in the thirties independently by Subranaman Chandrasekharcm and by Lev 
Davidovich LandauJj Their treatments differ enormously in style. 

After publishing his result ia an article, Chandrasekhar gave his deriva- 
tion in a 343 page classic bookEj which summarized: a) some of the results 
originally presented in German in the classic textbook Gaskugeln by Emden, 
b) basicAspects of quantum mechanics, developed in the classic book of R.H. 
FowlerEfl and in the work of P.A.M. Dirac, and 3) the detailed computations 
of the mass-radius relation for white dwarfs and their agreement with obser- 
vations with the clear statement of the existence of a critical mass equal to 
1.45M . 

The work of Lev Landau is totally different in its essential approach: 
Landau also uses all the expertise which could be gained from the classic 
work of Emden and quantum theory, at that time just formulated, and the 
corresponding relativistic corrections. Similarly he reaches the conclusion of 
the existence of a critical mass for a white dwarf, but now, everything, in a 
dramatic paper of no more then three pages. There by the stroke of crucial 
shortcuts in the physics of self-gravitating degenerate matter in relativistic 
regimes, Landau reaches that unexpected conclusion. Landau's paper is not 
only beautiful for the masterful interplay of physics and astronomy, but can 
also be considered the first paper introducing a new style of work in the 
scientific analysis of astronomical phenomena. 

Another non-negligible difference must be mentioned distinguishing the 
work of Landau and Chandrasekhar: Landau did not believe the possible 
existence in nature of the new concept of a critical mass he had just obtained, 
while Chandrasekhar never revealed such a critical attitude towards it. 
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At the time of presentation of his work in Cambridge, Chandra met with 
very strong resistance and criticism from Arthur Eddington, so much so that 
his paper published in the proceedings of the Royal Astronomical Society was 
preceded by an article of strong criticism by Eddington himself which clearly 
stated, with reference to Chandra's work, "... but its physical foundation does 
not inspire confidence, since it is a combination of relativistic mechanics with 
non-relativistic quantum theory." 

I was always interested in understanding some of the reasons which mo- 
tivated Eddington in his criticism. After all, at the time, Eddington was one 
of the greatest experts in the theory of generaLrelativity and was the author 
of fundamental textbooks both in astronomylij and general relativity:E3 for 
me more than his expressed criticism toward Chandrasekhar's work, it was 
interesting to understand whether there was any fundamental theoretical is- 
sue adopted by Chandrasekhar which could be disturbing to the eyes of a 
relativist. 

The answer to Eddington's criticism came finally in 1974 from Chan- 
drasekhar himself in his outstanding article "Why are the stars as they are" 
presented at the Varenna Summer school on the Physics and Astrophysics of 
Neutron Stars and Black Holes.EJ 

Nevertheless reconsidering recently the issue of the white dwarf critical 
mass, I have found a basic unproven assumption in Chandrasekhar's work. 
The distribution of the degenerate electrons is assumed to be uniform in the 
stars and unaffected by Coulomb interaction with the nuclei all the way to 
infinite densities. Indeed in Chandrasekhar's work the value of the critical 
mass is reached for an infinite central density of the star. On this issue 
Landau had made two important remarks: 

1. In 1938, following the earlier suggestions by Hund, LandauEl pointed 
out "it is well known that matter consists of nuclei and electrons. Never- 
theless it can be shown that in bodies of very large mass, this usual "elec- 
tronic" state of matter can become unstable. The reason for this lies in the 
fact that the "electronic" state of matter does not lead to extremely great 
densities, because at such densities electrons form a Fermi gas having an 
immense pressure. On the other hand, it is easy to see that matter can go 
into another state which is much more compressible — the state where all 
the nuclei and electrons have combined to form neutrons. ". Later Yakov 
Borisovich Zel'dovicha pointed out, quantitatively, that this "neutron 
condensation" for inverse beta decay starts to occur for hydrogen already 
at 10 7 ^3-. Then Harrison et s£5 pointed out that all configurations of 
white dwarfs after the onset of the inverse beta decay are indeed unstable 
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and the critical mass is reached not at an infinite density but at precisely 
the finite density marked by the onset of inverse beta decay. This result 
was confirmed W the analysis of pulsationaLmodcs of White Dwarfs by 
Chandrasekhar (see also Harrison et aE3). 

2. The second contribution of Landau is contained in the treatise with 
Yevgeny Lifshitz, where LandaulHI gives an estimate of the effects of the 
electromagnetic interaction between nuclei and electrons and explicitly 
shows how these interactions can be neglected only in the asymptotic 
regime for very large densities. The fiducial densities where indeed the 
Coulomb interactions can be neglected occur at p > 10 7 g/cm 3 for Z=56, 
namely at densities too high if compared to the density regime at which 
the stable white dwarfs and their critical mass occur. The interactions 
between nuclei and degenerate electrons therefore are not negligible for 
white dwarfs: they were certainly difficult to evaluate at the time of the 
Chandrasekahr's original work. 

Actually an interesting alternative approach, with respect to the evalua- 
tion of the interaction term between the electrons and the nuclei, was advanced 
two years earlier than the work of Chandrasekhar and Landau by the Soviet 
physicist Yacav Ilich Frenkel.Ej This paper, quoted by Chandrasekhar in his 
classic book,Ej has been in the past and is still today surprisingly ignored. 
Frenkel proposed to use a relativistic Thomas-Fermi model within a Wigner- 
Seitz approximation in order to describe stellar matter. The astrophysical 
motivations were not clear to Frenkel, though the theoretical formulation he 
proposed would have deserved a much more thorough examination. 

An important turning point came in this field due to the fundamental work 
of Feynman, Metropolis and Teller. While at los Alamos, they considered the 
equation of state of comnijessed matter.Q 

It was Ed Salpeteaffl .who first applied the Frenkel equations and the 
Feynman-Metropolis-Tellem approach to the treatment of white dwarfs: 
Salpeter's treatment was somewhat affected by the intricacies of the numeri- 
cal integrations and the drastic approximations he adopted in the numerical 
solutions. 

We ourselves rediscovered the Frenkel Equations and the relativistic 
Fermi- Thomas treatment applied to an appropriate Wigner-Seitz problem 
following the study of self-gravitating systems C3 and the riformulation of 
the treatment of these systems by using^ a_xelativistic generalization of grav- 
itational Fermi-Thomas like equations □ 113 lj. Today with relatively sim- 
ple numerical computations, we can reconsider the problem and present a 
comprehensive treatment by first considering the issue of the electromagnetic 
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interactions between—nuclei and electrons and turning later to the issue of 
inverse beta decaysJlil The new treatment£9 leads to the disappearance of 
all scaling laws with the chemical composition of the stars, which have been 
commonly assumed in all treatments up to now. A new value of the critical 
mass of white dwarfs is obtained, smaller than the one originally proposed 
by Chandrasekhar and by Salpeter, and one which strongly depends on the 
chemical composition of the star. These results appear in principle to be ob- 
servationally testable. Particularly attractive is the possibility of explaining 
the observed masses of binary pulsars: these systems are clearly neutron stars 
which could be formed by the onset of instability of stellar cores ven£_close 
to the critical mass of white dwarfs as presented by ChandrasekharO The 
evaluation of their masses will be clearly affected by our computations. 

3 The Relativistically Compressed Thomas-Fermi Atom 

In the well know.11 Feynman-Metropolis- Teller treatment of the classical Fermi- 
Thomas modelEnS an accurate numerical analysis is made of atomic configu- 
rations submitted to external pressure. Following Ferreirinho et al. u and 
Rufhni and Stella et o/.t3 we reconsider such a treatment in the framework of 
a relativistic generalization of the Fermi-Thomas model, which is here briefly 
derived and solved for different atomic species in different states of compres- 
sion. Results so obtained are applied to determine the equation of state of 
cold degenerate stars (white dwarfs) to show how the relativistic and elec- 
tromagnetic effects acting on microscopic scales affect the overall structure 
of the star. Finally, the equation of state obtained is used to determine the 
equilibrium configurations of white dwarfs. A comparison is made with the 
classic work of Chandrasekhar and Salpeter to show the differences introduced 
by the use of the new equation of state, and to show how those treatments 
correspond respectively to a zero-order and first-order approximation to the 
one we present here. 

Let us consider the spherically symmetric problem of a nucleus with Z 
protons and A nucleons interacting with a fully degenerate gas of Z electrons. 
The equation of electrostatics for this problem reads, with usual meanings for 
the symbols appearing in it, 

AV(r) = 47ren e (r) - Anen p (r) , (I) 

where we have introduced the number density of electrons n e (r) and of protons 
n p (r). 

In Eq. (|l|) the quantities V(r) and n(r) are of course not independent, 
being related by the equilibrium condition for a relativistic gas in a Coulomb 
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external potential (see 



0) 



c yP 2 F + m 2 c 2 — eV(r) = const = Ep , (2) 

where Ef stands for Fermi energy of the electrons. 

To put the Eq. (|l|) in a simple and dimensionless form we introduce the 
new function ^(r), related to the Coulomb potential by 

$(r) = V(r) + E F /e (3) 

and the corresponding dimensionless function x(r), implicitly defined by 

Zex 



$(r) = 



(4) 



Similarly we have the dimensionless expression for the radius r = bx, where 
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Let us express the second term of the source of Eq. (Q) in terms of dimen- 
sionless quantities: we assume a homogeneous nucleus, with a radius given by 
the approximate formula 

r nuc = 1.2A^ 3 Fermi (10) 

The number density of protons is therefore 

3Ze 



nuc 



9 (x r , 



(11) 
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Finally we can write Eq. (PP) in the form 



fx 

dx 2 



,3/2 
A/2 



1 



z 



4/3 



3/2 



■ix 



' ® (Xnuc Xj 



(12) 



where x nuc is the dimensionless size of the nucleus (r nuc = bx nuc ). Eq. ( |l2| ) 
is what we call the 11 generalised dimensionless Fermi- Thomas equation." 



The first boundary condition follows from x °^ r ^ an d therefore x 
and so 

x(o) - . 

The second condition comes from the normalization condition 
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with atom size ro = bx Q . Using Eq. (|l§) we obtain 
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which gives the simple relation 

N = Z[xox'(x )-x{xo) + l] • 
For a neutral atom we have N — Z and the Eq. (|l^) reads simply 

XoX'(x ) = x(xo) 



(13) 



(14) 



(15) 



(16) 



(17) 



Note that the physical quantities of interest, such as the Coulomb poten- 
tial and the density of electrons do not show any singularity in the center, nor 
on the surface of the nucleus, since they only depend on the function x an d 
his first derivative. The only discontinuity appears in the second derivative of 
X due to our rough assumption of a homogeneous nucleus. 

4 The equation of state 

Under these hypothesis, we determine the equation of state of compressed 
matter simply by computing for each value of the atomic compression param- 
eter xo the corresponding values of pressure and mass density: the first is 
generated by the degenerate gas of electrons and the second is given by the 
mass density of nuclei. We assume that the nuclei are arranged in a Wigner- 
Seitz lattice, each cell being filled by a relativistic gas of degenerate electrons. 
The shape of the lattice cell is unimportant in the lowest approximationJH so 
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we will assume a spherical cell, with a radius equal to the atomic radius. The 
average density of cells is thus 

_/ \ Am 3Am . . 

p{xo) = v^ = ' (18) 

where A and mo are respectively the number and mass of the nucleons. 

In this scheme the pressure is also easily determined, being simply the 
pressure generated by a Fermi gas of number density equal to the number 
density of electrons at the surface of the atom n e (xo), which can be expressed 
in terms of dimensionless quantities as follows 
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Our treatment differs from the zero-order approximation of Chan- 
drasekhar (uniformly distributed electrons) and also from the first-order ap- 
proximation of Salpeter in which the Thomas-Fermi equation for each value 
of Z and Xo was not solved and instead an expansion of the number density 
of electrons was introduced 



n(x) = no + e(x) . 



(20) 



This expansion applies only when the electrostatic interaction within nuclei 
and electrons is small with respect to the Fermi energy of electrons, i.e. in 
the limit of very high densities. 



5 Equilibrium Configurations for White Dwarfs 

We can now compute the equilibrium configurations for white dwarfs. The 
results are shown in figures ^ and ^|. We obtain, for the same central density, 
less massive configurations than those obtained by both Chandrasekhar and 
Salpeter. Namely for a fixed mass density, i.e. for a fixed value of A and of 
the compression parameter, our pressure is systematically lower than the one 
generated by a uniform shell of electrons. 

Numerical integration of the models is performed here for stars of oxygen 
and iron in the framework of a Newtonian theory of gravitation, as well as 
in General Relativity. We show that the first one is actually a fairly good 
approximation up to central densities of the order of « 10 9 g/cm 3 . 

It is shown that deviations from classical models are very important at 
the lower densities, up to 30% from the Chandrasekhar model and 10% from 
the Salpeter one. 
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Figure 1. Equilibrium configurations curve for Oxygen WD in the M/M@ — p c plane, ob- 
tained in Newtonian theory and General relativity. For comparison we phow the results 
obtained with Chandrasekhar model .Reproduced from Bertone & RuffmiM. 



We stop the integration at the critical density corresponding to the onset 
of inverse beta decays and we thus find the maximum mass for equilibrium 
configurations. This critical configuration marks also the onset of gravita- 
tional instability (see Harrison et al. EZl). 
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